N75-19590 


TUSKEGEE INSTITUTE 
SCHOOL OF MECHANICAL ENGINEERING 


TUSKEGGE INSTITUTE, ALABAMA 



CONTRACT NO. NAS9-13844 
BY 

JOHN P. SELLERS 
SPONSORED BY 

Reproduced by 

NATIONAL TECHNICAL 
INFORMATION SERVICE 

US Department of Commerce ■ 

Springfield, VA. 22151 | 

ENVIRONMENTAL AND THERMAL SYSTEMS SECTION 
MANNED SPACECRAFT CENTER 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
HOUSTON, TEXAS 

FECES -SUBJECT TO m'Z 

\ 


DECEMBER, 1974 




I 


NOTICE 

THIS DOCUMENT HAS BEEN REPRODUCED FROM THE 
BEST COPY FURNISHED US BY THE SPONSORING 
AGENCY. ALTHOUGH IT IS RECOGNIZED THAT CER- 
TAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RE- 
/ LEASED IN THE INTEREST OF MAKING AVAILABLE 
AS MUCH INFORMATION AS POSSIBLE. 



FOREWARD 


This report was prepared for Johnson Space Center of 
the National Aeronautics and Space Administration. The work 
was performed under Contract NAS 9-13344 with Mr. W. Ellis 
and Mr. B. French providing NASA guidance.* 

The work was performed from January 1974 to August 1974. 
The author would like to acknowledge the assistance of Mr. 

B7 Carroll, Lockheed Electronics Company, Inc. He would 
also like to express his appreciation to Tuskegee Institute 
for granting permission to undertake the task. 


ii 



ABSTRACT 


The analysis, begun in an earlier grant, of the 
NASA data obtained on a VCHP heat-pipe radiator system 
which was tested in a vacuum environment was continued. 
The study included further analyses and interpretation 
of the steady-state results and an initial analysis of 
some of the transient data. Particular emphasis was 
placed on quanitative comparisons of the experimental 
data with computer model simulations. Although the 
results of the study provide a better understanding of 
the system which was experimentally investigated, they 
did not, unfortunately, provide a complete explanation 
for the major test anomolies, namely, the observed low 
YCHP performance and the relatively flat radiator panel 
temperature distribution. The results of the study also 
suggest, for future implementation, hardware, software, 
and testing improvements. 
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1.0 INTRODUCTION 


As a result of L. B. Johnson Space Center sponsored 
feasibility heat-pipe radiator testing, valuable experimental 
data are being accumulated. In the study described herein, 
the analysis, interpretation and documentation of these data 
is continued. Particular emphasis has been placed on the use 
of computer thermal models. 

- An earlier study focused on the steady-state perfor- 

mance in the feasibility tests. Much of the data, however, 
were obtained under transient conditions. In the ultimate 
application of the systems envisioned, both transient and 
steady-state operating modes are important. Consequently, 
considerable emphasis was placed in the present study upon 
the transient results; and, in fact, distinguishing between 
steady-state and transient data. The latter effort was 
assisted greatly, but rather late in the contract period, 
by the computer program HPTRAN . 
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2.0 OPERATIONAL ANALYSIS OF FEASIBILITY YCHP 

The purpose of this portion of the study was to iso- 
late and examine phenomena pertinent to the operational 
characteristics of the feasibility VCHP. Particular attention 
is paid to its observed reduced maximum capacity. First, the 
complete feasibility VCHP results were examined, and certain 
data were selected, which seemed the most enlightening, for 
special study. 

The control equation for a VCHP with a cold-wicked 


reservoir is* 


= 1 + 


_r ( P _ P )/T _ 

V c V V /r v c 


< p v- p v> 


Since V R /V and m R /V are fixed design parameters , then 
from Eq . (2-1), it can be see# that the active VCHP condenser 


length parameter V' depends upon Py, T^ , and the inactive 
condenser length temperature T . It was shown in reference 
.1 that the effect of T on V' is small, so that T is largely 
determined by Py and T^ . 

The vapor temperature Ty corresponding to Py can be 
expressed bys 



T 

mean 


^REj / 01 


where T refers to the mean temperature of the coolant as 

mean 


it flows through the heat exchanger. Thus, l i' can either be 
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increased or decreased by changes in the inlet temperature 
T in of the coolant to the heat exchanger, the reservoir 
temperature Tr , and Cl. The conductance Cl depends upon 
the flow rate m of the coolant. 

The above brief summary is not limited to any specific 
wick and is deduced from a study of the f equation, Eq. (2-1); 
which is theoretically derived and is believed to be valid. 
Previous analysis (1,2) of the operation of the VCHP in the 
feasibility tests have shown that under design conditions the 
experimental ¥ values were less than theoretical when ¥ 
was about 0.4 or higher, resulting in a serious reduction in 
the heat transport of the VCHP. In fact, at maximum design 
conditions the actual ¥ reached a value only slightly over 
0.6, considerably less than the expected fully-open value of 
1.0. As a result, the effective area of the panel was only 
about one-half of its actual area. 


2 .1 Steady-State Evaluation 

The VCHP parameters that best define the operation of 
a VCHP are the operating temperature Ty, the heat rejected, 
Q rej , the reservoir temperature Tr , and the environment 
temperature T^, , (or Qp, Useful operating curves for a VCHP 
can, therefore, be obtained by various plots using these para- 
meters. Figures 2.1 and 2.2 taken from data used in reference 
1, are examples, where Ty is plotted versus and ¥ . 
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FIGURE 2.1 T y vs. Q REJ i % = T R and V R /V c =5-0 


VAPOR TEMPERATURE, T 


540 



490 * ‘ ^ -I II I1 1 I1TM— I I III fi l l 1^^^— ^.1 

0 .2 .4 -6 .8 1.0 

^a/l-c “ ^ 

FIGURE 2.2 T y vs. V j = T R and V R /V c = 5.0 
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Lines of constant Tg , (T^ assumed equal to T^ ) , are also 
included in the diagrams. The data presented are analytical 
values and were obtained from the steady-state heat-pipe 
radiator system computer program PRFORM ( 3 )* The principal 
characteristic of these curves is that for a given environ- 
ment the increase in operating or vapor temperature of the 
VCHP is small («8°) as the ammonia vapor - gas front moves 
from the VCHP fully-off to the VCHP fully-on positions. In 
addition to this small temperature rise, it can be observed 
that between the two limiting positions there is nearly a 
straight-line variation in vapor temperature. Therefore, 
a good indication of the operation, and control, of a 
VCHP during a test can be established from a plot of the 
VCHP vapor temperature versus the front location when compared 
with theoretical predictions. Unfortunately, in reality, the ■ 
front is not flat, and its position cannot be measured 
directly. In addition, the vapor temperature is not normally 
measured due to sensor installation complications. The front 
location. can be estimated from the axial outside surface temp- 
erature measurements along the VCHP condenser. The vapor 
temperature should be approximately the same as the surface 
temperature measurement of the "low - k" section located 
between the VCHP evaporator and the VCHP condenser. Using 
those data as the best available approximations of Ty and 
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<J ; , an attempt was made to construct operating curves, 

Figs 2-3 and 2.4, for the VCHP used in the feasibility tests. 

It was desired that each data point in a particular sequence 
would have nearly equal environmental temperatures and nearly 
equal reservoir temperatures, which, in addition to the 
necessity of steady-state heat transfer, greatly reduced the 
number of data points that could be used from the feasibility 
experimental program. Consequently, a ’'minimum'' of three 
data points could be found for only four test sequences. 

These were sequences 6, 7 (combined) and sequences 1 6, 17 
(combined). The former was made with a low design environment , 
and the latter with a high design environment. 

2.1.1 Sequence 6, 7 

Referring to Fig. 2.3 and Table 2 . 1 , as the inlet 
temperature to the heat exchanger was increased by a large 
increment, from 80 to 125 F, with T g = -83 to -90F, and 
T r = -83 to -86F, i.e., T R T £ , the vapor temperature, 

Fig. 2.4, and y changed remarkably little, 6 l to 62 F and 
0.43 to 0.45, respectively. Theoretically, the VCHP should 
have been fully open, *f' = 1 . 0 , at an inlet temperature 
slightly above 90F. The theoretical VCHP opening temperature, 

= 0, for sequence 6, 7 conditions is about 47F. It was 
pointed out in reference 1 that the feasibility VCHP actually 
opened a few degrees lower than predicted. Thus, it can be 
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Seq. No. 


Run Time 



6,7 

191-17-36 

13 

83 

6l 

.43 

6,7 

191-18-05 

15 

91 

61 

.43 

6,7 

191-18-45 

15 

97 

61 

WB 

6,7 

191-19-13 

17 

105 

62 

.45 

6,7 

191-19-33 

17 

125 

62 

• 45 

16,17 

192-07-51 

[ 30 

70 

70 

.32 

16,17 

< 192-10-50 

34 

72 

68 . 5 

.36 

16,17 

193.08-30 

54 

72 . 

71 

.42 

16,17 

192-23-25 

39 

77 

76 

• 56 

16,17 

192-07-31 

* 

30 

80 

79 

• 51 

16,17 

192-06-45 

30 

91 

79 

■ 58 

16,17 

192-07-06 

30 

94 

80 

. 6 

16,17 ' 

194 - 07-20 

73 

107 

86 

• 65 


TABLE 2.1 Feasibility VCHP data 


? 
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SEQ. 16&17 EXP. 



SEQ.6&7 EXP. 



SEQ, 6&7C AlC 




concluded that the cross-over point for the VCHP was between 
an inlet temperature less than 47F and the next experimental 
point at 8 OF and * 1 ' between 0 and 0.4. The cross-over point 
is where the experimental Ty became higher than theoretical. 
With an environment temperature of -83 to -90F, although 
the data are not complete, it can be deduced that the feasi- 
bility VCHP began to malfunction probably before it was l/4 
open. The large increase of Ty with T^ is typical of a 
wick that is in the process of drying-out, and in this case, 
prematurely. 

2.1.2 Sequence 16 , 17 

Referring to Fig. 2.3 and Table 2.1, as the inlet 
temperature to the heat exchanger was increased from 70 to 
72F, T e = -13 to -20F and T R = -13 to -19F, i.e. T R ~ T £ , the 
vapor temperature, Fig. 2.4, increased only slightly, and 
l V increased appreciably, from O .32 to 0.42, indicating 
good VCHP control. As was increased further to 107F, 
however, the vapor temperature rose sharply to a value of 
86 F v/ith «p not exceeding O. 63 . The VCHP opened at about 62 F. 
Thus, the cross-over point for the VCHP was between 62 and 
71F. With an environment of -13F to -20F, the data indicate 
that the feasibility VCHP was probably operating reasonably 
close to expectations up to an inlet temperature of 70F, 
or until the VCHP was about 2/5 open. This was some 


2-10 



improvement over sequence 6,7 but above 70°F again the 
large increase of Ty with T IN is typical of a wick that is 
drying out . 

2.1.3 Heat Transport Watts and Watt-inches 

In reference 1, the heat rejection capacity of the 
feasibility YCHP was represented by the parameter 
A better parameter for expressing the heat transport of a 
heat pipe is the parameter » where is the 

effective transport length and is determined for the feasi- 
bility VCHP from: 


L eff - * L o + * L ev 

was calculated from the feasibility steady-state 
data, and the highest values from the available data for the 
various sequences are presented in Table 2.2. Also presented 
are the values. It is interesting, and of possible 

significance, that Q jv^j approximately the same for each 
sequence , but Qp ^ j . L g ^ ^ increased as the environment temper- 
ature became larger . The implication here is that the VCHP 
watt-in transport capacity improved at the higher environments 
and operating temperatures. At the highest environment temper- 
ature, which was off-design, the VCHP performance apparently 
exceeded theoretical, due to higher than expected panel fin 
temperatures, (See Section 4.0.) At the maximum design 
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Zl-z 


SEQ. 

NO. I 

RUN TIME 

m 

lb/hr 

t in 

°F 

T v 

■°F 

6,7 

191-19-13 

290 

105 

1 q 

15 

192-15-30 

300 

106.5 

1 8 ^ 1 

16 

192-17-35 

H 33 

90 

it 

«. 

17 

193-01-10 

2090 

96 

87 

20 

192-12-00 



2000 

120 

114 

20 

193 - 11 - 15 ' 

1990 

109 

108 

21 

193-15-I8 1 

2000 

139 

139 


hr-ft 


REJ 

Btu 


hr 


1250 


1321 


Q REJ ,L eff 
Watt -in . 


* 

22 



.48 

1306 ! 

9.219 

58 

.62 

1293 

10,496 

59 

•63 

1256 

10,251 

60 

1 

. 66 

1311 

11 ,009 


11 ,563 


13,161 


TABLE 2.2 Q REJ and ^REJ'^eff Resu l‘ts for the Feasibility VCHP 













































































environment, sequence 1?, the VCHP maximum capacity was much 
less than theoretical. 

2 . 2 Transient Evaluation 

It is the intent now to point out, and investigate 
in depth, selected transient data that are believed to 
characterize the VCHP operation during the feasibility tests, 
and hopefully shed some additional light on its sub-performance. 
Equation (2-1), and the summary presented above, provide 
the theoretical background needed for the study. 

2.2.1 Tjpj Transients 

As previously discussed, and also in reference (l), 
the inlet temperature, Tj^, of the coolant at the heat 
exchanger entrance has a strong influence on Ty and conse- 
quently on V . Thus, the coprect measurement of T t , t is 
important if the data are to be compared with analytical 
predictions. In the feasibility tests the inlet temperature 
was measured with an immersion thermocouple installed just 
upstream- of the entrance to the heat exchanger. This measure- 
ment was designated AJ0020, Main Tube Inlet. In addition, 
a surface measurement, EN 001, was made on the inlet tube 
to the heat exchanger. Although AJ0020 should have provided 
the more accurate value of T-^ and has been used herein as 
such, there were periods during the test that EN 001 was 
apparently giving a better transient record of T^ N , if the 
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recorded data are valid. For example, refer to the experi- 
mental data for time period 192-06-06 to 192-06-14, Fig. 2.5 
and 193-03-00 to 193-P5-24. Figs. 2.6 and 2.7. In all of 
these periods , fluctuations of the VCHP temperatures can 
he seen which cannot he explained by reading AJ0020, since 
it remained fairly constant. The VCHP temperature traces, 
however, correlate very well with the recorded fluctuations 
of EN 001. If these data are valid, some doubt is cast on 
the use of AJ0020 for the coolant inlet temperature, parti- 
cularly under transient conditions. 

Early in the test series (191-15-08), the feasibility 
VCHP was inadvertently subjected to a near step change of 
of relatively large magnitude, Fig. 2.8. In about two 
minutes Tj^ was increased from a value of 55°F to l40°F and 
then after two minutes was decreased to 130°F. In addition, 
as the inlet temperature increased, the flow rate also went 
up from 200 lb/hr to 700 lb/hr for two minutes and then was 
returned to its initial value. The effect of the change 
in coolant flow rate is discussed below. Thus, the VCHP was 
subjected to an extreme thermal shock, which resulted in 
some interesting and unusual transient behavior. Because of 
the fast thermal response of the heat exchanger , when the 
inlet temperature was increased at 191-15-08 'the vapor temper- 
ature , sometimes called the VCHP operating temperature , rose 
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at nearly the same rate reaching a maximum value of 75 F and 
then dropping to 70F when the flow rate was decreased. The 
data seemed to indicate that if the flow had not been reduced 
would have reached a value of 0.6 in about A minutes. 

As it was, T achieved a value of 0.55 in about 3 minutes. It 
is interesting to note that of the six feeder heat pipes 
feeder C was most sensitive (indicating lowest capacitance?) 
to the sudden thermal load applied to the panel. The VCHP 
responded quite rapidly to the thermal shock, but it only 
opened about half-way, even though T IN was held between 125F 
to 150F for 11 minutes. This was certainly positive evidence 
early in the test program that the feasibility VCHP was not 
operating in the manner intended. 

Start-up of the feasibility VCHP by increasing T^^ 
was discussed in detail in reference 1 and also in Section 
3.0 of the present report. During shutdown the movement of 
the vapor-gas interface is toward the evaporator and the 
transients that occur are again of interest. An example of 
the shutdown of the VCHP by decreasing T-j-^ is presented in 
Fig. 2.9 where data are presented for the time period 
192-17-50 to 192-18-25. At the beginning the VCHP is open 
with V ~ 0.6, T jn - 100F, T r = -8F, = 60 Btu/hr - ft 2 

and m = 1200 lb/hr. During shutdown m increased to 1900 Ib/hr 
but this change had negligible effect on the shutdown. (See 
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the effect of m on W below). The inlet temperature was dropped 
from 100F to 75^ i- n about 5 minutes, (which resulted in 4* 
decreasing to slightly less than 0 . 5 ) held at 75 F for about 
5 minutes and then was decreased to 50 F in about 5 minutes, 
followed by a final decrease to 45F, after approximately 
another additional 18 minutes. The VCHP closed, X V = 0 , 
when Tj^ became 45F. 

2.2.2 T r Transients 

It was pointed out above that V can be made to in- 
crease or decrease by variation of , all other conditions 
remaining constant. When the feasibility VCHP was not at 
its maximum capacity the predicted effect of on H 7 was 
qualitatively confirmed by the experimental results. Referring 
to the experimental data, for the time period 193 - 07-43 to 
I 93 -O 8-361 Fig. 22-10, an elapsed time of 48 minutes, it can 
be seen that the inlet temperature, and rh are nearly 

constant at 71F and 1850 lb/hr., respectively. The reservoir 
temperature, however, in that period increased from - 26 F to 
-20F, or AT r ^ 6 F. The observed change in V was from 0.45 
to about 0.4 orA l F= 0.05* Although qualitatively correct, 
the observed value of A‘F/At^ was considerably less than 
expected. The vapor temperature, TC 3» remained nearly con- 
stant during this period, which was as expected. 
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Interesting VCHP operation occured between times 
I 9 O-O 9 - 3 O and 190-14-00, Fig. 2.11. In that period the 
VCHP was opened by the reservoir temperature decreasing 
while all other parameters were maintained constant. 

Setting l /' = 0 in Eq. (2-1), and from the experimental 
data, when T =0, Ty = T-j-^ = 65°F, T g — -50°F. For the feasi- 
bility VCHP, V R /V C = ?.5 and m R /V = 1.8. Using these 
data in Eq. (2-1), and solving for T R gives -6. OF, which 
compares very well with the experimental value of T R “* -7.0°F. 

After 190-14-00, the reservoir temperature continued 
to decrease, with and rii constant, dropping to a value of 
-80.0°F. At this point the observed W value was about 0 . 3 , 
much less than the corresponding calculated value of unity. 

2.2.3 m Transients 

Changes of the coolant flow rate may or may not affect 
the VCHP heat transport. For example, consider time 192-06-15, 
Fig. 2.5- The flow rate of the coolant was decreased from 
a value of 1950 lb/hr, the approximate design flow rate, to 
1100 Ib/hr at 192-06-18. During that change the reservoir 
temperature increased slightly, which would tend to decrease 
V , but Tjn was constant at 84°F. The major perturbation 
was the 50 per cent decrease in flow and it did not result 
in any decrease of Ty from its initial value of 78°F or any 
detectable drop in the VCHP condenser temperatures. The 
latter implies that a 50 per cent reduction in the design 
coolant flow did not affect l !' . 
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On the other hand, the data of Fig. 2.8, which were 
used in the discussion above to illustrate the transients 
during a large, rapid increase of also provide an illus- 
tration where a decrease in A did affect the VCHP operation. 

At 191-15-10 the flow was lowered from 700 Ib/hr to 200 lb/hr, 
which almost simultaneously lowered the vapor temperature 
from 75F to 73F, and Y from over 0*5 to between 0.4 and 0.45* 

The data of 192-22-47, Fig. 2.12, illustrate the 
thermal response of the VCHP when the flow rate is rapidly 
and drastically reduced. From an initial value of 1840 Ib/hr, 
rn was reduced to below 400 lb/hr, and held at that level for 
two minutes. The flow was then momentarily shut-off before 
being brought back to its original value of 1840 Ib/hr. 

This entire flow excursion took a little over ten minutes. 
During this period a significant drop in the readings from 
VCHP thermocouples nos. 5-13, that is, those surface thermo- 
couples on the active portion of the VCHP before the flow 
was reduced, can be seen. 

Consider now the data for time period 192-16-00 to 
192-18-40, Fig. 2.9. During this time period the reservoir 
temperature was -8F and constant, and the VCHP was opened 
twice by increasing The only apparent difference between 

the two was the magnitude of the coolant flow to the heat 
exchanger: 1150 Ib/hr and 1875 lb/hr for the first and second 
start-ups, respectively. The lower flow-rate start-up began 
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at a vapor temperature (equal to the inlet temperature) 
perhaps as much as 7° lower than the start-up with the higher 
flow rate. It should be noted that in this sequence, as the 
inlet temperature was increased the VCHP condenser tempera- 
tures were decreasing at the point in time when the VCHP 
turned on. 
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3.0 LOCKHEEDS TRANSIENT HEAT PIPE HPTRAN PROGRAM 


3 • 1 Need for a Transient Program and its Use in the 
Present Contract 

Regardless of the particular type of radiator being 
considered for heat rejection, its ability to thermally res- 
pond within certain time limits to a change of operating 
conditions is an important characteristic of the system. It 
is not unreasonable to assume that during many duty cycles 
the heat rejection apparatus will be operating more frequently 
in a transient mode than in steady state. 

A purported advantage to be gained by utilizing heat 
pipes in a heat-rejection system is a relative fast thermal 
response, but in order for the designer to make meaningful 
comparisons of candidate systems, the transient character- 
istics of each must be available and studied. 

The feasibility heat-pipe radiator utilized a cold- 
wicked reservoir which, compared with wickless reservoirs, 
has a faster theoretical response, but again the lack of good 
quantitative data for each for comparison is evident. 

There is a need for both experimental and analytical 
transient heat-pipe data, particularly under simulated 
operating conditions and obtained from experiments incorpor- 
ating full-scale apparatus. The results obtained from the. 
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feasibility heat-pipe system can be used as an excellent 
starting place toward fulfilling this particular need, 
specifically, for the case of a cold-wicked reservoir. 

Although much data were taken over a wide range of 
operating conditions during the 102.3 hr. feasibility test 
program, the test variables were changed frequently through- 
out the test and it was often difficult to establish' those 
times when steady-state had actually been achieved. Conse- 
quently, only about 10 steady-state points, Table 3*1» were 
listed in reference 1, an initial analysis of the test 
results made by the present author. With the aid of the 
transient computer program HPTRAN , (5)» an objective under 
the present contract was to review and search the same 
data for additional steady-state points which, if successful, 
would permit a broader and more complete analysis of the 
results than was carried out in reference 1. Under the 
present contract this objective was not completely achieved 
for two reasons: 

a, The transient computer program, HPTRAN, was not 
completed until late in the contract because of revisions. 
This greatly reduced the amount of available time for its 
use. 

b. Programming difficulty was experienced converting 
the program to the computers available in the Tuskegee area. 
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3 • 2 Checkout of HPTRAN 


Soon after the initiation of the present contract, 
Lockheed Electronic Corporation, Houston, Texas, at the 
request of NASA, began work on a transient heat-pipe radiator 
systems program. The program, which they named HPTRAN, (5)> 
was an extension of their earlier steady-state program PRFORM. 

Before HPTRAN could be used with confidence in the 
present investigation, however, it was necessary to make 
certain checks and to establish that the program gave rea- 
sonably valid results over the range of conditions of interest. 
The experimental data ( 4 ) obtained in the feasibility tests 
could be used in the comparisons and validation. 

Most of the initial comparisons between HPTRAN results 
and the experimental data gave very poor agreement. It was 
not certain how much of the disagreement was due to the new 
transient additions to this program and how much was due to 
the fact that the position of the vapor-gas front in the VCHP 
condenser of the feasibility hardware did not agree with the 
analytical predictions (1,2). The latter discrepancy was 
known prior to the writing of HPTRAN and .is due to the sub- 
performance of the VCHP. 

It was clear that the only method of checking the 
transient portion of HPTRAN using the feasibility test data 
v/ as to input into the program the experimental front position 
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versus time. Unfortunately, the actual front is not flat 
and thin as assumed in the program and its location as well 
as width cannot he determined directly from the experimental 
data. They can,however, he estimated from the VCHP experi- 
mental surface temperature measurements. Using those data, 
the approximate front location as a function of time was 
established and, after a program modification, could be 
input ed to HPTRAN rather than calculated from the program 
equations. With the front then positioned fairly well in the 
program, a comparison of the computed and experimental panel 
temperatures provide a reasonable check of the transient 
operation of the heat pipe system, but primarily for the 
feeder heat pipes and the panel . Since the actual front 
location at a given time may itself depend on the transient 
heat transfer within the VCHP, this approach of inputting 
to the program the experimental front position does not 
necessarily provide a check of the computed VCHP transients. 

The validation of HPTRAN was the subject of reference 
6. Carroll, the author, who also wrote the program, concluded 
that HPTRAN does simulate transient response of the feasi- 
bility VCHP radiator system fairly accurately. In his 
study, emphasis was placed on the transients associated with 
starting-up of the VCHP from an M off" condition by increasing 
the temperature of the coolant entering the system. 
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Table 3*2 lists the particular time sequences which 


Carroll used in his program validation. Actually an additional 
time sequence was included but was omitted in Table 3-2 since 
the reservoir temperature was not correct. 

The present investigator agrees with the conclusion 
of reference 6 , that the program HPTHAN does indeed have 
merit, and since this was substantially shown in reference 6 , 
the arguments need not be repeated here. Whether the program 
has been satisfactorily validated is questionable, in the 
opinion of the present investigator. In that regard it would 
be worthwhile to take the runs of Table 3.2, for example, 
and investigate each of them in greater depth than was done 
in reference 6 . 

3-3 Results for Time Period 193-07-03 to 193-08-36 

To illustrate, consider the time period 193-07-03 
to 193~08-36 of Table 3*2, a span of 92 minutes, where Tj^ 
was increased from 57°P to 72°F by one large initial rise 
and later an additional small increase, Fig. 3*1* The reser- 
voir temperature also went up from - 39 F to -20F, but at a 
near constant rate, Fig. 3*1* The heat absorbed, , and the 
coolant flow rate were nearly constant at 56 Btu/hr-ft^ and 
18?0 lb/hr. , respectively. This was a sequence 17 run, and 
in fact, includes the data at time I 93 -O 8-30 which were used 
as a steady-state point of Table 3*1 ■ 
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TABLE 3 . 2 Time 


FLOW RATE 
Ib/hr 


DURATION 

MIN. 



3iC 


■> 


1 58-80 

-70 

285 

64 

| 52-66 

-56 

1924 

41 

J 63-72 

-16 

10?0 

47 

| 57-72 

-27 

1870 

92 


es Used for Program Validation 
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The calculated front location expressed by the dimen- 
sionless parameter, V' ca j c » and the front location as 
determined from the VCHP condenser surface temperature mea- 
surements, */' , are also presented in Fig. 3.1. These 

are the V' values used in reference 6 for the time period 
193-07-03 to 193-08-36. Referring to the ^ calc curve > 
which is the analytical front location as calculated from 


HPTRAN , according to reference 6 it is calculated from; 


f (^-V P v" 

W c / — Tj 


-Pw ) 


m g R g 


where 


! Np 

£ (P V“ P V S )/T s 

N = Number of nodes on the inactive portion of 
p the VCHP 


i = Node number on the inactive portion of the VCHP 
. In the earlier steady-state program PRFORM, (3), 

Y was obtained from a somewhat different equation; 



where T g is the average temperature of the Inactive portion 
of the VCHP. As the VCHP opens during the time period 
under study the V calc curve, Fig. 3.1, rises to a value 
of 0.82 due to the increase of the inlet temperature (in 
spite oP the VCHP closing effect of an increasing reservoir 
temperature), and then both the T^ N and Y . calc curves have 
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leveled-off at the 455 minutes mark. Henceforth, the y , 

/ calc 

curve does not appear to he correct on the basis of the 
and data of Fig. 3*1* For example, from time 440 to 475 
Tj^ is constant, but increases from -29°F to -25°F. On 
the basis of Eq . (2-1), and presumably Eq . (3-1), should 
decrease accordingly after 440 minutes but with some time 
lag due to heat storage in the transient mode. A significant 
effect of on ¥ was found in the steady-state study of 
reference 1 for nearly the same conditions as in the present 
discussion. See Fig. 3.2. It would appear that in the 
HPTRAN calculation of ¥ calc reserv °i r temperature 

effect has been attenuated, since f c alc does n °t decrease 
in the 35 minute period but remains constant. It is diffi- 
cult to explain this result on the. basis of the transient 
operation and further study is needed. 

If Qpgj rad is the net heat radiated from the panel, 
then a heat balance for the entire system gives: 

Q IN = Q REJ = ^Stored + Q out = Q Stored + Q REJ rad. 3-2 


30 

0UT )=C 2 (T i, 

i=1 


In Eq . (3-2) the numbers are nodal point designations. At 
steady state the heat storage term is zero. Thus, 
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28 


C p A ^ T IN 


T ) = 
1 OUT 



REJ 


- Q 


rad 


REJ 


In the transient mode Eq. (3-1) may he rewritten: 


Q stored ~ C p m ^ T IN T OUp QREJ ra d QrEJ QREJ rad 

It can he seen in Fig. 3*1 that at time 510 

minutes (measured from 193-00-00) calculated by HPTRAN 
using ¥ e has leveled-off at a steady-state value of 
813 Btu/hr. The Q^-gj value which was calculated from the 
experimental panel temperatures listed in reference 1, 
at this time was 980 Btu/hr., and in reasonable agreement 
with HPTRAN. 


Another check of the numerical programs can be 
made at the 510 minutes point. Referring to Fig. 3*1* 
the value based on ^ ca q c from HPTRAN is about 1320 

Btu/hr., and y ca ]_ c = 0.84. The analytical value given in 
reference 1 and calculated by the steady-state program 
PRFORM (3) is 1238 Btu/hr at ¥ calc ( = 0.?4. Ostensibly, 
the differences between the two programs can be explained 
by: 


a. Some differences in the input data to each of 
the programs, particularly, T ^ , T^ and m g R g/ v c • 

b. Some improvements in the analysis contained 
in HPTRAN but not PRFORM. 
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c« The ¥ calc from HPTRAN being suspect as dis- 
cussed in the preceding section. 

Consider now the transient operation during the 
start up of the VCHP and referring once again to Fig. 3 •! » 
it can be seen that the Q REJ curves are characterized ini- 
tially by severe fluctuations and a relative large peak 
value, 1200 Btu/hr for Q REJ ( Y EX p). As is °P ened 

by a rather rapid increase of Tj^ , one would indeed expect 
a surge of heat flow into the cold system and a large initial 
value of Q rej , but the fluctuations are believed to be due 
to the numerical method used in the program, and could likely 
be reduced by increasing the number of nodal points inside 
the VCHP. This should be considered, as a greater number 
of nodal points may also significantly improve the accuracy 
of the results in the transient mode. 


Figure 3.3 compares the calculated VCHP and experi- 
mental axial wall temperature distribution for three 
different times, initially (423) » 28 minutes later (451 )» 
and at the end of the designated time period (516). The 
calculated values are the ¥ results from HPTRAN. It 

can be seen that at time 516 minutes the experimental 
temperatures were lower than the calculated values beyond 
28 in. along the VCHP. The experimental data indicate 
that the VCHP temperatures between 20 inches and 34 inches 
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VCHP TEMPERATURE, 



FIGURE 3.3 VCHP Temperature vs. Distance for Run No. 193-07-O3 
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decreased after 464 minutes. This decrease which was 
not predicted by HPTRAN and may be due to the movement 
of the front toward the evaporator, heat conduction, or 
heat leaks, is also brought out in a plot of the VCHP 
temperatures versus time, Fig. 3*4. The measured decrease 
of the VCHP temperature at 24 inches (TC 14) is about 15°F 
in 50 minutes where the HPTRAN results gave only a 3 degree 
drop! 

The panel temperatures versus time are presented in 
Figs. 3-5 a nd 3 . 6 . It can be seen that the root temperature 
of Feeder D, which begins to drop after 480 minutes, is 
the only measured panel temperature indication of the 
T r closing effect which was quite evident in the VCHP data 
and discussed in detail above. The panel temperature 
distribution at the end of the time period is shown in 
Fig. 3.7 where the peak points on the temperature distri- 
bution curve are the root temperature of the active 
and partially active feeders. The "valley” minimum points 
are located half-way between feeders. Section 4.0 of the 
report is a study devoted to finding an explanation of 
why the temperature difference between these two locations, 
i.e., on the feeder and midway between, is less experi- 
mentally than calculated. 
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PANEL TEMPERATURE, 



TIME, min. 

FIGURE 3.5 Panel Temperature vs. Time for Run No. I93-07-03 



PANEL TEMPERATURE, 







4.0 PANEL TEMPERATURE DISTRIBUTION 

Comparisons made in references 1 and 2 and Section 
3.3 of this report "between the experimental results obtained 
from the feasibility tests and the computer data have brought 
out that a consistent disagreement exists in the temperature 
distribution on the panel fin between active feeders. That 
is, considering two adjacent active feeders, the difference 
between the two high panel temperatures, located on the 
f ee g er heat pipes, and the minimum temperature at the mid- 
point position between feeders , was consistently less than 
the calculated temperature difference for corresponding points 
on the panel . 

The comparisons are presented for the runs of Table 

3.1 in Figs. 4.1 - 4.6 and Table 4.1, where the computer 
data were obtained from the steady-state program PRFORM. 
Referring to Table 4.1, it can be seen that the theoretical 
temperature difference T nR F j^p T F^ T p is several degrees greater 
than experimental. For example, by 3 to 6 degrees for feeder 
A, by 6 to 13 degrees for feeder B, and by 5 to 9 degrees 

for feeder C. The magnitude of A T p depends upon the quantity 
of heat Q being conducted away from the feeder heat pipe 
through the panel. Thus, to insure that Q 0xp ^ Q calc . 
consider run 193-15-18 where Q REJ gxp^REj ca ic . by 
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TEMPERATURE, 



193 - 14-50 



DISTANCE, in. 


FIGURE 4.‘l 


Panel Edge-to-Edge Temperature Distribution* 
Time 193-14-50 ’ 
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TEMPERATURE , 



FIGURE 4.2 Panel Edge-to-Edge Temperature Distribution; 

Time 191-18-4.5 
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T E MPERATURE, 



■ FIGURE 4.3 Panel Edge-to-Edge Temperature Distribution- 

Time 191-21-40 




temperatures 


19 2-15-30 

0 


0 

0 EXP. 

0 

A CALC. 

193-08-30 


DISTANCE, in. 

FIGURE 4.4 Panel Edge-to-Edge Temperature Distribution; 

Times 192-15-30 and I 93 -O 8-30 
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193 - 10 - 20 
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o EXP. 
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DISTANCE, in. 

FIGURE 4.6 Panel Edge~to-Edge Temperature Distribution; 

Times 193-10-20 and 193-15-18 
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NA = Feeder not fully active. 
TABLE 4.1 Panel Fin-Length /\T . 
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20 percent (see Table 3«1)« From Table 4.1 for run 193-15-18, 
AT? FHp ) A T p FHp by 3° for feeder A, by 9° For 
feeder B and by 7° for feeder C. It can be concluded, there- 
fore, that the measured panel AT between the feeder and the 

x'' 

midpoint between feeders was significantly less than pre- 
dicted by the analysis. 

The radiation from a point on the panel varies with 


the fourth power of the local temperature so a few degrees 
of panel . temperature inaccuracy can have a large effect on 
Qp EJ . Since the panel temperature distribution is crucial 
to the proposed optimization study of the panels, an effort 
was made under the present contract to understand these 
experimental results and to make any corrections to the 
analysis found necessary. 

The afore-mentioned discrepancy between the analyti- 
cal and experimental results was surprising, since the prob- 
lem is unrelated to the efficiency of operation of heat pipes 
and would appear* to be simple heat conduction in the panel 
with radiation to and from the panel surface facing the 
simulator. The opposing surface of the panel was insulated 
in the tests and is considered an adiabatic in the computer 
program PRFORM . 

PRFORM results are obtained from a model consisting 
of a nodal point located on each of the feeders plus three 
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nodal points equally spaced between feeders. One-dimensional 
heat transfer perpendicular to the feeders was assumed. For 
convenience, in the present study of the panel temperature 
distribution a separate small computer program was written 
to exepedite a parametric study of an isolated section of the 
panel located between two active feeders. In this small 
program, titled SNOD, originally the number of nodal points 
between feeders was set at seven, four more than in PRFORM. 

In addition, results obtained from a computer program FIN 1 (7) 
(9 nodal points) were compared with the experimental data. 

4. 1 Effect of Heat Transport in Feeder, Q 

Figure 4.7 shows that AT^ increases linearly with Q. 
Since Q for an individual feeder cannot be measured experi- 
mentally, it is difficult to compare analytical and experi- 
mental results for the same Q value. For this reason, it 

is desirable that Q„ YT3 > to conclude that AT < 

HjAl caj_c . p exp 

AT p calc.' This was illustrated above for run 193-15-18. 

4. 2 Comparison of Results from Different Programs 

The panel temperature distribution calculated from 
three computer programs, PRFORM, SNOD and FIN each using 
different numerical methods are compared in Fig. 4.8. The 
conditions chosen correspond to run time 193-10-20. It is 
seen that good agreement exists between the results from 
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q' a = 93Stu/hr-ft 2 

k = 95 Btu/hr -ft -^F 
t = .02 in 


'TIME 1 93-10-2 0 


O EXP, 



.7 


.8 


.9 


• -3 .4 .5 .6 

NOND1MENSIONAL FIN DISTANCE, x/L f 

FIGURE 4.8 Comparison of Panel Temperatures Calculated from j Different Computer Programs 


PRFORM and SNOD, but FIN produced a flatter temperature pro- 
file which gives A ca p c i n closer agreement with the 
experimental results, but before any conclusions can be 
drawn Q must also be compared. The feeder heat-transport Q 
is input in SNOD but is calculated from the output from the 
other two programs : 

Q - C 2n (Ty - T m ) PRFORM 

Q = -kA ( — ) . FIN 

dx x = 0 

Q calculated from the FIN temperature curve, Fig. 

4.8 is unrealistically low, 112 Btu/hr compared to 1?0 Btu/hr 
for the other two programs. It would appear, therefore, that 
the actual temperature distribution is approximately as de- 
picted by the sketched- in dashed curve of Fig. 4.8. Unfor- 
tunately, an insufficient number of temperature measurements 
were taken to draw the actual temperature distribution. 

4.8 Effect of Panel Thermal Conductance on At 

— p calc . 

An increase of the panel thermal conductivity from 
85 to 100 Btu/hr-ft-F (18 percent) reduced AT^ ca p c from 
8 to 7°F or 1°, Fig. 4.9* The nominal value of k is 95 
Btu/hr-ft-F. It is reasonable to assume that the true value 
of k is within 1 Qfo of this nominal value. Thus, it is con- 
cluded that an incorrect value of k in the programs cannot 
explain their disagreement with the experimental results. 
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An increase of the panel thickness from 0.0185" to 
0.0215’' (16 percent) reduced A ca jc from 7.8 to 6. 9°F 
or about 1°F, Fig. 4.9. A nominal value of 0 . 020 ” is used 
in the calculations. 

4.4 Effect of QJ on AT 

• — A p calc . 

An increase of with Q constant increases the mag- 
nitude of the panel temperatures but has little effect on 
A Tp calc ^ , Fig. 4.10. In the normal program calculation it 
is assumed that is uniform across the panel. An investi- 
gation was made to determine the effect on AT n of a 

p calc . 

higher value midway between feeders than at the feeders. 
For simplicity, a straight line variation of was assumed. 
The results are shown in Fig. 4.10. For example, a variation 
of from 82-5 to 92.5 Btu/hr-ft 2 gave little reduction in 
ATp cal0i compared to when was constant at 93 or 80 
Btu/hr-ft 2 . 

4.5 Effect of Panel Insulation 

The panel simulator was wrapped in insulation after 
its installation in the large vacuum chamber. A simplified 
analytical attempt was made to determine if the insulation 
could explain the observed temperature distribution. 

Referring to Fig. 4.11, 

^ ~ ^ simulator + ^ insulation 

- A h rs (T n - T e ) + Ah ri (T n - CT e ). 

- where C = Ti/Tg 
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PANEL TEMPERATURE, 


540 


k = 95 Btu/hr f t F 
t = .020 i n 
Lf = 4 in 

SNOD PROGRAM 


530 - 


Q a = 82. 5-92. 5 Btu/hr ft 2 


(4T p =12°) 


( AT p =12.6°) 


(AT p - 1 3.6°) 


0 1. 2. 3. 4. 

DISTANCE FROM ACTIVE FEEDER, IN 


FIGURE 4.10 Effect of Absorbed Heat Flux QJ on Panel Temperatures 



At low environment temperatures, in most runs, it is 
expected that the insulation temperature exceeded the envir- 
onment temperature, hence C>1.0, perhaps 1.0 - 2 . 0 . At 

high environment temperatures, T. « T„ and C =» 1.0. 

1 Ji 

It can be seen from Fig. 4..12 that AT decreases 

P C3.XC « 

as C increases. Unfortunately, the flat temperature distri- 
bution observed in the tests occurred at high environments 
as well as low environments. Thus, the insulation does not 
appear to be the explanation for the observed panel temper- 
ature profile. 

4 • 6 Effect of Fin Length 

If the radius of the feeder heat pipes is considered, 
the distance from the feeder outside wall to the midpoint 
between feeders-, or fin length, is reduced from 4 " to 3 . 69 ", 
Fig. 4.13. The shorter fin length when used in the FIN 
program reduced A T p by about 12$. Thus, some of the dis- 
crepancy can be eliminated by considering the true feeder 
envelope rather than an imaginary line heat source as is 
presently assumed in the computer models. 

a_ 

4 . 7 Effect of Panel € Ratio 

Using the FIN program a /€ was varied from 0 to 1.1 
( 6 = 0 . 9 ), and the results are presented in Fig. 4.14. It 
can be seen that A T p calc # decreases as a/e increases . 

Normally this ratio is taken as unity. It is unlikely that 
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FIGURE 4.13: Feeder Pipes and Fin 



a/e 

FIGURE 4.14 Effect of Panel a/e Ratio on Panel Fin - Length At 
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a/€ was greater than unity in the feasibility tests. Thus, 
the experimental panel temperature distribution cannot be 
explained by e^a* 

4 . 8 E ffect of Number of Nodal Points 

The number of nodal points between the feeder and 
the midpoint between feeders was increased from three to 
nine. It is seen Table 4.2, that the increase had negligible 
effect on AT 

p calc . 
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Number of 

Distance 

from Active Feeder, in. j 

Nodal Points 

o 5 

2 

4 

Five 

. 

540R 

532 . 1R 

529. 6r 

j Eleven 

j 540R 

j 532. 2R 

529. 6r 

' 


TABLE 4.2 Comparison of Panel Temperatures 
for Five and Eleven Node Models 








5-0 CO NCLUSIONS 


a. With an environment of -83° to -90°F, the fea- 
sibility VCHP began to malfunction before it was 1/4 open. 

b. With an environment of -13° to -20°F, the data 
indicate that the feasibility VCHP was operating satisfac- 
torily up to inlet temperatures of 70°F or until the VCHP 
was about .4 open. 

c. The feasibility VCHP watt-in. transport capacity 
improved at the higher environments and operating temperatures 
At the highest environment temperature, 105°F, which was 
off-design, the VCHP performance exceeded the analytical 
value . 

d. The inlet temperature measurement, AJ 0020, is 
suspect, particularly for transient conditions. 

e. ...An estimate of the maximum vapor-gas front 
movement in the direction of increasing the active portion 
of the VCHP condenser as a result of a severe thermal shock 
was from 0 to 3/5 fully-open in about 4 min. The indications 
were that the thermal response of a correctly functioning 
VCHP would be relatively fast as expected. 

f. Chronologically, the VCHP was malfunctioning 
from the outset of the test program. 
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g. The VCHP was Q RE j limited, but not Q REJ • f 
limited. This could be considered evidence that the VCHP 

low heat transport problem was centered in the heat- exchanger- 
evaporator rather than in the VCHP condenser. Additional 
evidence was that at the design coolant flow rate the VCHP 
vapor temperature exceeded the coolant heat exchanger out- 
let temperature. 

h. At this stage of its development the HPTRAN 
computer program has merit and is useful, but requires 
further validation. Some of this can be accomplished using 
the existing 'feasibility test data. For example, the H' 

calc. 

results for time- period 193“ 0 ?-°3 to 193-08-36 are suspect. 
Final acceptance of the program awaits experimental data 
obtained on a VCHP with gas-vapor front locations that 
agreed much better with theoretical predictions. 

i. - The. difference between the maximum and minimum 
panel temperatures between two active feeder heat pipes 
was consistently less in the feasibility tests than calcu- 
lated. A detailed parametric study of the panel temperature 
distribution did not disclose any one reason for this test 
anomoly . 

,J. The computer models should be modified to con- 
sider the true feeder envelope and thereby a shorter panel 
fin length value than is currently used. 
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6 . 0 REC OMMENDATIONS FOR FUTURE STUDY 


Based on the results of the studies described in 
this report, the following recommendations for future study 
can be made: 


1. Determine why the feasibility YCHP was 

limited but not limited. 

2. The experimental transient behavior discussed 
in depth in Section 2.2 should be compared with the com- 
puter results from HPTRAN . 

3. Continue with the validation of HPTRAN in the 
manner illustrated in Section 3 » 3 « When the need for 
greater computer accuracy is justifiable, the number of 
YCHP nodal points should be increased. 

4. Continue with the panel temperature distribution 
study to try to' bring the analytical and experimental 
results into closer agreement. 
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8 . 0 SYMBOLS 


A 

C 

C 

c 


p 

Cl 

C2 

h 


eff 


ca 


ev 


m 



Q 


a 

Q 


A 

REJ 


Area 

Capacitance 

V T E. 

Specific heat of the coolant 
Conductance 

Conductance ' 

Radiation heat transfer coefficient 

Thermal conductivity of panel 

La. + Lev , 

2 2 

Length of the active portion of the VCHP 
condenser 

Length of the VCHP evaporator 

flow rate of the coolant 

mass of the non-condensible gas 

Pressure of the vapor in the VCHP 

Partial pressure of the vapor in the VCHP 
reservoir 

Partial pressure of the vapor in the inactive 
portion of the VCHP 

Heat transport in a fully active feeder heat 
pipe 

Absorbed heat flux 

Total .heat removed from the coolant 
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Q 

< 3 . 

R 

{ 

t 

T 

■T 


rad 


g 


n 


nR 


T. 

l 


-T 


IN 


T 


OUT 


T 


mean 


T, 


T 


E 


T 


T 

i R 

VCHP 

V R 


Net heat radiated from panel 

Heat transfer 
Gas constant 

Thickness of radiator panels 
Temperature of the panel at a nodal point 

Temperature of the panel at a nodal point 

located on an active feeder heat pipe 

Insulation temperature 

Temperature of the coolant as it enters the 
heat exchanger 

Temperature of the coolant as it leaves the 
heat exchanger 

Logarithmic mean temperature bet ween the coolant 
in the heat exchanger and Ty 

Temperature of the vapor in the YCHP 

Temperature of the environment 

Temperature of the vapor and gas in the inactive 
portion of the VCHP condenser 

Temperature of the reservoir 

Variable conductance heat pipe 
Volume of the reservoir 


V c Volume of the vapor in the VCHP condenser 

x Distance across the panel perpendicular to the 

feeder heat pipes 

CL .. Absorptivity of the radiator panels 
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T nR T rimin 

Emissivity of the radiator panel 
Stefan-Boltzman constant 

Ratio of the active length to the total length 
of the VCHP condenser 
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